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Abstract — Frequency Estimation of a complex exponential is 
a problem relevant to a large number of fields. In this paper 
a computationally efficient and accurate frequency estimator 
is presented using the guaranteed stable Sliding DFT which 
gives stability as well as computational efficiency. The estimator 
approaches Jacobsens estimator and Candans estimator for 
large N with an extra correction term multiplied to it for the 
stabilization of the sliding DFT. Simulation results show that the 
performance of the proposed estimator were found to be better 
than Jacobsens estimator and Candans estimator. 

Index terms : Frequency Estimation, Sliding DFT, Radar 
Signal Processing 

I. Introduction 

The frequency estimation of a complex sinusoid is a fun- 
damental signal processing problem which finds a lot of 
importance in fields like radar signal processing Jn radar signal 
processing the computational capability not only needs to be 
accurate but also fast which involves estimating the frequency 
every second. The received signal is corrupted by white 
gaussian noise. The best coarse frequency estimation of the 
signal is from the peak of the N point DFT of the received 
signal. A N point DFT is typically calculated for a data 
length of N samples which gives a resolution of ^ . For real 
time spectral analysis a well known computationally efficient 
method is the sliding DFT especially in the cases when a new 
DFT spectrum is needed every few samples.The sliding DFT 
is computationally efficient than the radix-2 FFT. The sliding 
DFT performs a N point DFT within a sliding window of N 
samples.The window is then shifted by a sample for the next 
iteration and a new N point DFT is calculated which utilises 
the old N point DFT values. The single bin SDFT algorithm 
is implemented as an IIR filter, to which a comb filter can 
be added so as to compute all N DFT spectral components.In 
practical applications the algorithm can be initialised with zero 
input and zero output.The output wont be valid until N samples 
have been processed. The sliding DFT though has a marginally 
stable transfer function because all its poles lie on the z- 
domain's unit circle.The filter coefficient numerical rounding 
might force the poles outside of the unit circle which results 
in instablity.A damping factor r may be used to force the 
pole into a radius of r inside the unit circle which guarantees 
stability.But the value of the bins obtained from the guaranteed 
stable sliding DFT are different from that of the traditional 
N point DFT. So the main problem lies in a fine estimation 
of frequency from the guaranteed stable sliding DFT which 



differ from the traditional N point DFT. The fine frequency 
estimation generally follows the coarse frequency estimation 
which is in turn done from the N point DFT.But the resolution 
of this process depends on the spacing of points taken for the 
fine resolution.In [3] -[6] the fine resolution estimate is done 
through a function on the DFT bins estimated from coarse 
resolution frequency estimaton.In [7] a simple relation for DFT 
interpolation is used and in [8] a bias correction is provided 
for [7]. In this paper we derive Jacobsen's formula for the 
guaranteed stable sliding DFT for which a bias correction term 
is also derived which includes bias correction for the damping 
factor r as well. 

II. PROBLEM DESCRIPTION 

A single complex sinusoid with white gaussian noise can 
be represented in the form 

r[n] = Ae^'^'' ^ w[n] 



where A and uj are unknown variables which represent the 
amplitude and frequency of the complex sinusoid respectively 
where ^nd kp is the index of the peak of the 

sliding DFT. 6 is to be estimated from the three samples around 
the peak of the sliding DFT where \5\ < 1/2 

The transfer function for N point sliding DFT filter can be 
represented as 
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where |r| < 1 Therefore a particular output bin of sliding 
DFT is written as 
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where |r| < 1 In practice though a particular output bin 
can be found out using the following recursive relation which 
basically serves the computational efficiency purpose 

Rk[n] =Rk[n- l]re^'2^^/^ - x{n - N)r^ + x{n) 

The sliding DFT bin where the peak occurs and its immediate 
neighbours can be represented as follows:- Let the indices for 
the peak be kp and that of its inmiediate neighbours be kp — 1 
and kp-\-l respectively. 
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R[kp] = AJ2 r"e^T^''" + w[kp\ 



n=0 



Af{5) + w[kp 
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R[kp - 1] = ^ ^ ^n^j^{S+l)n ^ ^[^r.^ _ 
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= Af{S + 1) + w[kp - 1] 



(2) 
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R[kp + 1] = A ^ j."eJTf ("5-1)" + w[kp + 1] 
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where w[A;] is the DFT of w[n] which also is white and 
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f{a) = r"e^-^' 



n=0 

Our aim is to estimate the value of S from these three 
samples R[kp],R[kp — 1] and R[kp + 1] so that 

cj = 27r/N{kp^6) 

becomes the fine frquency estimate. The two stage process 
consists of finding kp in the first stage and S in the second 
stage. 

III. PROPOSED ESTIMATOR 

To determine S from the set of three equations we consider 
the Geometric Progression sums of each of the DFT bin and 
solve for S using the approximation that the second and higher 
powers of 5 are negligible as compared to S. We basically 
exploit the Geometric Progression sum in this case because of 
the added damping factor r which in turn satisfies the relation 
|r| < 1 and thus makes the common ratio in the Geometric 
Progression to be less than one. 

AT-l 



1 — re^ N 



where \r\ < 1 To estimate S We evaluate the first differ- 
ence and second difference of f{S) For evaluting the first 
and second differences which are f{S-\-l) — f{S — 1) and 
f{S + 1) - 2f{S) + f{S - 1) respectively (1) is used. Using 
elementary trigonometry the first difference 

/(,5 + l)-/(,5-l) 
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which can be further simplified to 

2jrsin(27r/Ar)e^'2^'^/^(l - r^) 
- 1 + ^2gi47r5/Ar 2rei2^VA^ cos(27r/Ar) 

Using elementary trigonometry the second difference 

f{5 + l)-2f{5) + f{5-l) 
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.27r(5 + l) ^ T 27rS ' .27r(5-l) 

1 V I - re^ n 

- (1 + ^2gi47r5/iv 2re^'2^^/^ cos(27r/7V))(l re^'27r5/^) 

We get the following relation for the ratio of the two 
differences 

/((5+l)-2/(^) + /((5-l) 
_ -jcot(7r/Ar)(l - - 2jr s\ii{2t:5 /N)) 



l + r2 + 2rcos(27rVA^) 

For large N and 6'^ <<1 we use the following approximations 

sin(27r^/Ar) = 2ti5/N and cos{2tt5/N) = 1 the above 
relationship can be simplified to 

_ -j cot(7r/A/')(l - - AjrnS/N) 

_ 47Tr6 cot (tt/N) j ( 1 - r ^ ) cot (tt /N) 
~ 7V(l + r)2 (l + r)2 

Simplifying the above relationship further we get 

f{5 + l)-f{S-l) ^ 
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f{s^i)-2m^f{s-iy 

AirrS cot(7r/A^) 



Ar(l + r)2 

(4) Thus an estimate of S can be produced by the substitution 
of f{S) R[kp]J{S - 1) R[kp + 1] and f{S + 1) 
R[kp — 1] (which were inturn derived in the relations (1),(2) 
and (3)) in (4); 

. ^ + tM^[/Nl R[kp - 1] - R[kp + 1] 
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'2i?[A:p] - R[kp - 1] - i?[A:p + 1]' 



IV. Numerical Comparisons 

In this section we present a numerical comparison between 
the proposed estimator and the other estimators namely, Can- 
dace's estimator and Jacobsen's estimator.The performance of 
the proposed estimator is compared across various damping 
factors which is shown in Fig. 1. The simulation is done 
in MATLAB.A sinusoidal signal is taken whose frequency is 
varied from 30.1 MHz to 30.9 MHz in steps of 0.1 MHz. 128 
samples of the signal are taken where the sampling frequency 
is 128 MHz. A 128 point sliding DFT is taken for the proposed 
estimator while a 128 point traditional DFT is taken for the 
other estimators. The noise taken is white gaussian noise. 
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Average Peak Location Estimates (bin) vs Bin I^umber 




Average Peal< Location Estimates (bin) vs Bin Number 
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Fig. 1. Performance of the proposed estimator with different damping factors 



Fig. 3. Performance of the estimators in presence of noise 



Fig. 2 shows the bias values of the proposed estimator and 
the other estimators with respect to the ideal frequency in 
the absence of noise for N =128.This shows that none of 
the estimators are unbiased.For the proposed estimator the 
damping factor taken is r =0.9 for the simulations.lt can be 
clearly seen that the bias of the proposed estimator is less than 
1000 Hz while those of the other estimators have higher bias 
values. 
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Fig. 2. Variation of bias of different estimators in absence of noise 

Fig. 3 shows the performance of the proposed estimator and 
the other estimators with respect to ideal frequency is presence 
of noise for A/^ =128. 

Fig. 4 examines the bias of the estimators in presence of 
noise. It can be seen that the proposed estimator has lower bias 
as compared to the other estimators. With higher SNR values 
the bias values in presence of noise approach the bias values 
in absence of noise. 

Fig. 5 examines the RMSE values of the estimators for large 
number of observations. 

The SNR values taken for the above simulations range 
from 2 dB to 3dB. At higher SNR values the estimator bias 



Estimator Bias vs Bin Number 




Fig. 4. Variation of bias of different estimators in presence of noise 

and the variance values decrease further and all the three 
estimators behave nearly the same with the proposed estimator 
giving even lower bias and variance values. It should be noted 
that with the damping factor r approaching 1 with higher 
values of N the proposed estimator have nearly the same 
performance.The value of the damping factor has been kept 
around 0.9 for guaranteed stability and better performance.lt is 
intutively satisfying that the bias correction factor in Candace's 
estimator is also a part of the derived expression for the 
proposed estimator.The fewer number of operations involved 
for the proposed estimator compared to that of other estimators 
and the lower bias and variance values makes the proposed 
estimator very useful for radar signal processing. 

V. Conclusions 

A new estimator is proposed which requires very few 
number of operations per output sample.The estimator has a 
correction term for bias. The good performance of the estimator 
is justified in the paper.The proposed estimator has low bias 
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Estimator Variance vs Bin Number 
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Fig. 5. Variation of RMSE of different estimators in presence of noise 

and variance values which makes it a really valuable tool in 
the field of radar signal processing. Simulation results showing 
superiority of the proposed estimator is provided. 

The present work revolves around accurately estimating a 
single tone frequency or a multi tone frequencies where the 
tones have large seperations.A potential future work is the 
extension of present work to accurately estimating multi tone 
frequencies which are closely seperated without increasing the 
number of computations. 
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